Instability of antiferromagnetic magnons in strong fields 
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We predict that spin-waves in an ordered quantum antiferromagnet (AFM) in a strong magnetic 
field become unstable with respect to spontaneous two-magnon decays. At zero temperature, the 
instability occurs between the threshold field H* and the saturation field H c . As an example, we 
investigate the high-field dynamics of a Heisenberg antiferromagnet on a square lattice and show 
that the single-magnon branch of the spectrum disappears in the most part of the Brillouin zone. 
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PACS numbers: 75.10.Jm, 75.10.-b, 75.50.Ee 

There are several reasons of studying the effect of 
strong magnetic field on quantum AFMs. A growing 
family of weakly interacting spin systems, which includes 
chain jTj, ladder Q, and square-lattice AFMs S, al- 
lows now experiments in previously unreachable field 
regimes. The high-field physics is proven to be rich for 
one-dimensional (ID) AFMs, to mention only incommen- 
surate gapless modes in a spin-i chain 0. At the same 
time, in 2D or 3D the field evolution of the Neel ordered 
ground state is trivial. Spins cant gradually from the 
antiparallel structure until the magnetization saturates 
at the critical field H c . Zero-point fluctuations vanish at 
the same field making the ground state at H > H c purely 
classical. One can suggest that similar quasi-classical sce- 
nario applies also for the spin dynamics. In contrast, we 
predict in this Letter that on the way to the saturated 
phase the excitation spectrum of an ordered AFM under- 
goes unexpectedly strong transformations. In high fields 
below H c , magnons are overdamped and disappear in the 
most part of the Brillouin zone (BZ). 

The effect originates from the field-induced hybridiza- 
tion of single-magnon states with two-magnon contin- 
uum. Previously, Osano et al. Q investigated the effect 
of such an interaction on the two-magnon dynamical re- 
sponse at low fields H <C H c . In this limit the hybridiza- 
tion leads to a weak renormalization of the one-magnon 
spectrum. In a strong field, however, two-magnon con- 
tinuum overlaps with the single-magnon branch and pro- 
duces instability of the latter. At T = 0, there is a thresh- 
old field H* (H* 0.76H C for square and cubic lattices) 
above which magnons become unstable with respect to 
spontaneous decays. The argument for that is based on 
softening of the spin-wave velocity of the sound mode 
near the AFM wave- vector Q. For small k = k Q, the 
magnon energy is approximately 



uj\s_ w ck(l + ak 2 ). 



(1) 



In zero field both the classical and renormalized disper- 
sions bend downward implying a < P-R]. On the 
other hand, at H = H c the low-energy asymptote for 
Uk coincides with the corresponding expression for a fer- 



romagnet: Uk « P. Simple continuity arguments predict 
that the high-field upward curvature of the spectrum for 
k — > changes to the low-field downward bend at some 
field H* < H c , where the cubic factor a vanishes. There- 
fore, a is positive in the field region H* < H < H c and 
long-wave magnons are kinematically unstable towards 
the spontaneous decay ||. 

The above argument for the existence of the thresh- 
old field H* is general and valid for all ordered quan- 
tum AFMs irrespectively of their dimensionality, length 
of spin, or position of Q. In the following, we study 
in detail the high-field dynamics of a spin-i Heisenberg 
antiferromagnet with nearest-neighbor interaction on a 
square lattice described by the Hamiltonian 



(2) 



(id) 



The ground state of this model is ordered at < H < H c 
for arbitrary value of the on-site spin S . 
For the canted AFM phase 



S? 



Sfe iQ - r i cos 6 + Sf sin 9. 
-_ S r < - - 



5? sin 6> 



,x e iQ ' r > cost 



(3) 



Sf" = Sf , and Q = the boson representation of 

H is obtained by applying the Dyson-Maleev transfor- 
mation in the twisted frame: Sf — S — a]^, S^ = 
V2S(1 - a\aj2S)a t , Sf = y/2Sa\. The tilting angle 6 is 
determined from vanishing of the linear boson term Tl^: 
sin 6 = H/8S, H c = 8S. Apart from a constant, the bo- 
son Hamiltonian is a sum of quadratic, cubic and quartic 
terms: H = H {2) + H (3) + H [A) Q. Neglected terms 
containing products of five and six boson operators are 
of higher orders in an expansion parameter cos 2 9/2zS ||, 
z being the number of nearest-neighbors, and give small 
corrections even for S = h- 



The quadratic Hamiltonian 
#00 = ^[Atoiok - ^k(a k a-k 



+ 4«- k ) 



(4) 



where A k = 45(1 + sin 2 6»7 k ), B k = 4S'cos 2 6'7k, and 
7k = I (cos k x -\~cos fcy), is di&gonalized by the Bogoliubov 



i 



transformation = Uk&k+w k 6^ k with u^+v^ = Ak/W, 
2itkWk = Bk/ujk, giving the classical spin-wave energies 



w k = 45v/(l + 7k)(l-cos20 7 k). 



contributions to the dispersion arise from 7v 4 ) and from 
the renormalization of the tilting angle Gof : 



(5) E (3) (k) = 4(u k + vl) [-n cos 29 + A cos 2 6> - m sin 2 6> 



The magnon spectrum is defined in the paramagnetic 
BZ. It has a sound mode near k = Q with the classical 
spin-wave velocity c = 2^/2S cos 9. 

In the leading order in 1/zS quantum corrections to 
the magnon dispersion can be found from the Dyson 
equation for the normal Green's function: G(k, t) = X) 
— i(T6 k (£)6 k ), neglecting anomalous contributions to the 
self-energy Q. To find E(k, uj) we express the interac- 
tion terms via quasiparticle operators &k and consider 
the first-order perturbation from the quartic part and 
the second-order perturbation from the cubic part. 



+7k(— m cos 29 + |<5 cos 2 6> — n sin 2 #)] — 8itk^k 
x [ \ Asin 2 9 - \ mcos 2 6» + 7 k ( Acos 26> - /kms 2 6» + ± Jsin 2 6>| , 
EW(k) = 8sin 2 0(A-n+m)[( Uk +« 2 )(l-7 k )-2 UkUk7 k]. 



where various two-boson contractions are defined as n = 
k vl, m = £k u k7k> 8 = Ek u k«k, A = Ek«k«k7k. 
Four-magnon vertices do not contribute to the dispersion 
to this order in 1/zS. 

Perturbation theory gives the magnon energy with the 
first 1/5* corrections as 



pert 



cck + 5> W (k,^k). 



(9) 





(a) 



(b) 



FIG. 1. Lowest-order diagrams contributing to the magnon 
self-energy. 

There arc two interaction terms with cubic vertices 
fromW^) 

^ = 5T E rg 3 (4 & t &1 + h.c.), 



l-2-3=Q 



(6) 



l+2+3=Q 



where we abbreviated momenta with their subscripts. 
The first term V\ is responsible for the hybridization be- 
tween single- and two-magnon spectra and is the princi- 
pal interaction of the problem. The vertices are given by 



F 



(i) 

1:23 



8Ssin20$(k; q), 



$(1;23) = 7i (ui +v 1 )(u 2 v 3 + v 2 u 3 ) +72(^2 + v 2 )(uiw 3 
+ V1V3) + 73(^3 + V 3 )(u 1 U 2 + U1U2), (7) 
(2) 1 — 

and r^23 = — v8S'sin2(9F(k, q), where similar expression 
for F(k, q) can be found in po[ . The lowest-order con- 
tributions to the self-energy shown in Fig. 1 are 



E (1) (k,w) =4S*sin 2 26>^ 



$(k;q) 



UJ - CJ q - Wk-q+Q + lO 



(8a) 



zW(k,uj) = -4Ssin 2 29Y . ( 8b ) 

The first term describes a virtual decay of a magnon into 
two-particle intermediate states. Frequency-independent 



Fig. 2 shows classical and renormalized spectra for 5 = \ 
in two representative fields. At H = 0, the only non- 
vanishing correction E^ 3 ^(k) yields a 16% enhancement 
of magnon energies in the entire Brillouin zone |^J^] . In 
low fields the magnon spectrum resembles zero-field re- 
sult except the vicinity of the BZ center. The mode at 
k = describes a uniform precession of spins about H 
and, for axially symmetric systems its classical frequency 
luq = H is exact, i.e., not changed by quantum effects 
Jl2| . With increasing field the cubic terms grow and push 
down magnon energies. Above 0.7 H c , correction from 
E^ (k, Wk) exceeds 50% of Wk signifying a break down of 
the 1/S perturbation expansion for spin- wave frequencies 
(^) . The correct renormalized spectrum plotted in Fig. 2 
is recovered by solving the Dyson equation 



Wk — E(k, uj) = 0. 



(10) 



Standard l/S expansion breaks down for any finite S 
when the bottom of two-magnon continuum becomes 
nearly degenerate with a part of the single-magnon spec- 
trum for H close to the decay threshold field H* . Un- 
renormalized "classical" value of H* is calculated by ex- 
panding Eq. (^|) up to the third order in k. The cubic 
factor a is anisotropic and vanishes first for k x — k y at 
H* = -^H c « 0.756ff c . 

At H > H * the magnon self-energy acquires imaginary 
part due to spontaneous two-magnon decays, which obey 
the energy conservation law 



Wk - Wq + Wk-q+Q- 



(11) 

Spin- wave is stable, if Eq. ([Ll]) has only trivial solutions: 
q = k, Q. In this case, the two-magnon density of states 



p 2 (k,w) = £ 5(g 



Wk-q+cj) vanishes at w < uj^. 



P2 (k, uj) has Van Hove singularities determined by sym- 
metry extrema k — q* + Q = q*+G, where G is a 
reciprocal lattice vector. At H > H* , one of these points 
q* = (k + Q)/2 crosses the one-magnon state if 



2 



uj k > 2w( k+ Q)/ 2 . 



(12) 



When applied to the phonon mode, Eq. (|1J) is equivalent 
to the previously used condition a > Ojnote the differ- 
ence between k and k in Eqs. (||) and At H > H*, 
classical magnons (||) decay in the region around k = Q 
and remain stable in the vicinity of uniformly precessing 
mode k = with the decay threshold boundary given by 
equality sign in Eq. (|12|). 
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FIG. 2. Magnon dispersion for k = n(r),r)): (i) thin dashed 
lines represent the classical spectrum o>k; (ii) dot-dashed lines 
are u£ cr ; (iii) solid lines are solution of the Dyson equation. 
For H = 0.1H C , curves (ii) and (iii) are indistinguishable. 

Pitaevskii argued that such boundaries are actually 
termination points of the spectrum ||. This fact is a 
consequence of logarithmic divergence of the one-loop 
diagram near the threshold: EW(k ( u) oc ln(-AQ/R), 
where Co = Aw + v c Afc, v c is a velocity of created 
magnons, and R is a cut-off. Termination point for the 
decay threshold into a pair of rotons is known in the spec- 
trum of superfluid 4 He. Since the self-energy is singular, 
a proper renormalization of the cubic vertex is necessary 
in that case fl. Our problem, however, is different, be- 
cause two new magnons created in an elementary decay 
process are again unstable. Near the decay threshold, 
the dominant contribution to E^ 1 ' arises from the close 
vicinity of q = (k + Q)/2, therefore we can approximate 
Imw q w Imwk_ q+ Q = |T. Even small T's remove com- 
pletely divergence of the diagram (Fig. la): 



SW(k,w)cx 



In 



\J~Klu 



7,2 _ 



r 2 



R 



7T Auj 

— + arctan — ^- 



(13) 



This analysis leads us to two conclusions: (i) vertex cor- 
rections can be neglected, but instead a self-consistent 
treatment of intermediate magnons and their decaying 
rates is required; (ii) the boundary ( |l2"| ) between stable 
and decay regions in BZ is smeared and at H > H* all 
magnons simultaneously acquire finite life-times. 



In the limit k — > the decaying rate is small and can 
be found analytically. The self-consistent Born approxi- 
mation (SCBA) with dressed lines in Fig. la gives 

ImS (1) (k,w) = -47rSsin 2 20 
d 2 q 



(2nf 



$ 2 (k;q)</>(w-w q -w k _ q+q ), (14) 



where function <p(u>) is normalized by J <f>(uj)duj — 1 and 
has a characteristic width of the order of a sum of decay 
rates of created magnons. Using the asymptotic form 



$(k;q) m -\[kq(k - ^/V^cos 3 ^ 1 / 2 we find 
w k w ck(l + ak 2 - z/3fc 2 ), 



(15) 



with = (9 sin 2 6/64^^/2 cos 3 6») 2 / 3 . Similar to a 3D case 
H , this expression for p is valid beyond a narrow field re- 
gion near H* determined by higher-order nonlinearities. 

Another region in BZ, where damping is small and ex- 
citations are well defined, lies near the zone center. This 
follows from the same arguments Jl2| , which predict no 
effect on luq from quantum fluctuations, since the exact 
frequency at k = is real and has no imaginary part. 

The spin-wave damping in the entire BZ is consid- 
ered by using the magnon spectral function A(k, uj) = 
— -ImGQz, oS), where G(k, uS) is calculated in the Born 
approximation. Results for A(k, oS) at H — 0.85H C and 
H = Q.9H C are presented in Fig. 3(a) and 3(b), respec- 
tively. In the non-self-consistent approximation (dashed 
lines), A(k, ui) consists of a narrow one- magnon peak and 
two-magnon side-band, which exhibits loose maximum at 
higher energies. The quasiparticle peaks survive even in 
the classical decay region Eq. ( |l2| ) , because hybridization 
pushes them out from the two-magnon continuum. This 
spurious feature arises due to the lack of self-consistency, 
as the two sides of Eq. ( |il"l) are treated with different ac- 
curacy in the non-SCBA. That is, the magnon energy is 
renormalized while the energy of the continuum (right- 
hand side) is given by a classical expression. Peaks disap- 
pear in the SCBA, which takes into account modification 
of two-particle continuum due to renormalization of one- 
magnon states. We performed numerical solution of the 
Dyson equation in the SCBA with 96 x 96 k-points in BZ. 
Results are plotted in Fig. 3 by solid lines. Our analy- 
sis shows that a self-supporting instability intensifies at 
H = 0.85H C for magnons from the region k < Q/2, which 
are classically stable but become strongly damped since 
they lie in the decay region for renormalized spectrum. 

Magnons peaks in j4(k, uj) are strongly suppressed 
(H = 0.85H C ) and disappear (H = 0.9H C ) because of 
the shift of single-particle pole in G(k, uj) from the real 
axis into the complex plane due to rapidly growing de- 
cay surface (|TT|). This is in contrast with Pitaevskii's 
scenario for superfluid He, where the pole's residue be- 
comes vanishingly small Single- magnon excitations 
reappear again only at H « 0.997J C , where the decaying 
vertex rW(l; 23) is small. 
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FIG. 3. Magnon spectral function for several points 
k = 7r(?7,7?), (a) H = 0.85H C , (b) H = 0.9H C . Thin dashed 
lines represent the non-SCBA. Solid lines are the results of 
the SCBA. Corresponding values of rj axe shown near each 
curve. 

To complete our study, we investigate now the high- 
field dynamical response of the square-lattice AFM, 
which is probed in neutron diffraction experiments. Our 
calculation of the dynamical structure factor S a P(k, ui) = 
J dt(S^{t)S^ k )e i,vt is standard. We express it at 
T = via the spin Green's function G a ^(k, t) = 
-i(TS^(t)S^_ k ) by S a/3 (k,uj) = -26(uj)G al3 (k, «). The 
latter is related to the boson Green's function by using 
Eq. (||) and the Dyson-Maleev transformation. The re- 
sult for the inelastic part of the dynamical structure fac- 
tor is 

S xx (k, u) = ttS sm 2 9(u u + v k ) 2 A(k, u) + it cos 2 (9 



X ^(ttqtfc-q+Q+WqUk-q+Q^Cw-Wq-Ufc-q+Q), 

q 

SW(k, u>) = nS(u k - v u ) 2 A(k, lu), (16) 
S zz (k, w) = ttS cos 2 (9(w k _ Q + Vk-Q) 2 A(k - Q, w ) 

+7T Shl 2 ^(u q tI k -q + UqMk-q) 2 (5(w -CJ q -CJk-q), 

q 

where we use the unperturbed spectrum for the two- 
magnon contributions. These two-magnon contributions 
described by the last terms in S xx and S zz are negligible 
in high fields because of small factors Uk ~ cos 2 #. In the 
same approximation Uk ~ 1 except for close vicinity of 
Q. Thus, in high fields H* < H < H c the longitudinal 
component S zz (k, ui) reduces to its elastic part, whereas 
transverse components S xx (k, ui) and S vv (k,uj) are pro- 
portional to AQc, ui) with approximately same momen- 
tum independent prefactors. 

In conclusion, the nonlinear coupling between one- and 
two-magnon states, which exists only in the canted AFM 
phase, becomes very important in high fields. Together 
with the field-induced kinematic instability of the spec- 
trum this interaction leads to suppression and disappear- 
ance of single- magnon excitations for H* < H < H c . 
Therefore, an intriguing situation arises in the high field 
regime: the ground state of an ordered quantum AFM 
is nearly classical, while the excitation spectrum and the 
dynamical response are strikingly different from the clas- 
sical results. Though, our analysis was restricted mainly 
to the spin-i AFM in 2D, the qualitative picture remains 
valid for other values of S and for 3D. 
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